Circularly polarised X-rays as a probe of non-collinear magnetic order in multiferroic 

TbMnOg 



F. Fabrizi,i'2 H. C. 
A. T. Boothroyd,^ 



Walker, ^'Q L. Paolasini/ F. de Bergevin,-^ 
D. Prabhakaran,'^ and D. F. McMorrow^ 



o 
o 

(N 
>>■ 



I 

o 
o 



> 

00 
O 
rn 

vn 
o 
o\ 
o 



European Synchrotron Radiation Facility, Boite Postale 220, 38043 Grenoble, France 
^London Centre for Nanotechnology and Department of Physics and Astronomy, 
University College London, Cower Street, London WCIE 6BT, UK 
^ Department of Physics, Clarendon Laboratory, University of Oxford, UK 

(Dated: May 19, 2009) 

Non-resonant X-ray magnetic scattering has been used to study the magnetic structure of mul- 
tiferroic TbMnOa in its ferroelectric phase. Circularly polarized X-rays were combined with a full 
polarization analysis of the scattered beam to reveal important new information on the magnetic 
structure of this canonical multiferroic. An applied electric field is shown to create a magnetic nearly 
mono-domain state in which the cylcoidal order on the Mn sublattice rotates either clockwise or 
counter-clockwise depending on the sign of the field. It is demonstrated how this technique provides 
sensitivity to the absolute sense of rotation of the Mn moments, and to components of the ordering 
on the Tb sublattice and phase shifts that earlier neutron diffraction experiments could not resolve. 

PACS numbers: 75.25.-l-z, 75.50.Ee, 61.05.cp 



Multiferroic materials exhibit unusual physical proper- 
ties as a result of coupling between the various forms of 
spontaneous ferroic order they display [JJ] . Recently, con- 
siderable interest has been generated following the dis- 
covery of a large magneto-electric coupling in TbMnOa 
, which was subsequently shown to be due to the estab- 
lishment of non-collinear antiferromagnetic order driving 
the formation of a ferroelectric state (FE) [3]. A cen- 
tral theme of ensuing studies in TbMnOa and other sys- 
tems has been to explore the electric field control of mag- 
netism, and correspondingly the magnetic field control of 
FE 0, S] ■ Of paramount importance in the development 
of any theory of the magneto-electric coupling in this 
class of materials, is a complete and accurate microscopic 
description of the magnetic structures they display. 

Historically the utility of circularly polarized X-rays is 
well established in dichroism experiments on ferromag- 
nets d, 0) Hi ■ Less attention has been paid to the case of 
antiferromagnets, with even fewer examples of their use 
in diffraction from non-collinear systems such as helices 
and cycloids [13] ■ Here we report on a new applica- 
tion of circularly polarized X-rays in combination with 
full polarimetry of the scattered beam to the refinement 
of magnetic structures in magneto-electric multiferroics. 
Besides the classical advantages provided by X-ray mag- 
netic diffraction (discrimination between spin and orbital 
contribution, in non resonant scattering) [111 ], we demon- 
strate how this technique provides unique insight into the 
formation of cycloidal domains in TbMnOs. A key fea- 
ture of these experiments was the control of the popula- 
tion of magnetic domains by an applied electric field. 

Bulk measurements first established that TbMnOa un- 
dergoes a transition below T—27 K to a multiferroic state 
that is both magnetic and FE, and demonstrated how a 
magnetic field could be used to control ferroelectricity 
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FIG. 1: a) Crystallographic and b) magnetic structures of 
TbMnOs. The incommensurate magnetic structure propa- 
gates along the b direction with periodicity 1/r foo ~ 3.55 60 
on both the Mn and Tb sublattices. Both the transverse and 
the previously undetermined longitudinal components of the 
Tb moments are shown, c) Projection in the b — c plane of 
the two cycloidal magnetic domains, also showing the newly 
determined longitudinal Tb moment component. Domain 1 
(2) is favoured by electric field E > (E < 0) (see Fig. E}. 



0. A major breakthrough in understanding this effect 
was provided by a neutron diffraction study which con- 
cluded that at T=27 K the Mn 3d magnetic moments un- 
dergo a transition from a coUinear to a non-collinear, cy- 
cloidal phase, described by an incommensurate wavevec- 
tor km — rb*, which removes a centre of inversion so 
as to form a FE state For this class of multifer- 
roic, the FE moment P may be conveniently written as 
P oc km X C, where C —^.^ Si x S.^+i characterises the 
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magnetic structure adopted by the spins Si H, 13, Ml- 
In the case of TbMnOa the Mn magnetic moments in 
the cycloidal phase rotate in the b — c plane (Fig. [T]), 
and with km parallel to the b axis, P is either parallel 
or anti-parallel to the c axis depending on the sign of C. 
The same study [1] proposed that the Tb 4/ moments are 
also sinusoidally modulated at k,„, but are transversely 
polarized along the a axis. Polarized neutron diffrac- 
tion has recently been applied to multiferroics [3, 
including TbMnOa [r?|, where it was found that the cy- 
cloid could be switched from propagating clockwise or 
anticlockwise, depending on the direction of the applied 
electric field 17]. Our non-resonant X-ray magnetic scat- 
tering (NRXMS) study of TbMnOa estabhshes the utility 
of X-rays for studying the electric field control of mag- 
netism in multiferroics. Although NRXMS experiments 
are in many ways much more demanding than neutron 
diffraction, due to the extreme weakness (~ 10^^ — 10^^) 
of the NRXMS relative to the charge scattering, they are 
nontheless of considerable value due to the complemen- 
tary nature of the information they provide. In the case 
of TbMnOa this has allowed us to reveal the ordering of 
a b component on the Tb sublattice and to determine 
various phase relationships that were hitherto unknown. 

TbMnOa adopts the Pbnm space group at room tem- 
perature with ao = 5.302 A, bg = 5.856 A and cq = 7.400 
A [13]. We selected a single crystal grown by the fioat- 
ing zone method, as used in previous resonant X-ray 
scattering experiments [Toj . The sample had dimensions 
2 X 2 X 0.8 mm'^, and a poling electric field E of about 
1 kV/mm was applied during the cooling procedure in- 
side a ^He evaporation cryostat. Two copper electrodes 
were glued by high conductive silver paste on the sample 
c-surfaces, and the electric field was removed during the 
X-ray experiments. 

X-ray magnetic diffraction experiments were per- 
formed at the ID20 beamhne 20] (ESRF, Grenoble, 
France) using a monochromatic beam at a wavelength 
of A =1.66 A, so as to excite non-resonant processes, 
which allows the spin and orbital magnetizations to be 
determined independently [2l[. (Neutron scattering, by 
contrast, is sensitive to their sum.) Linearlly polarized X- 
rays were converted into left-circular ly polarized (LCP) 
or right-circularly polarized (RCP) photons using a dia- 
mond quarter-wave plate (circular polarization typically 
99%), or allowed to pass unchanged through the phase 
plate, see Fig. [2ja). For a cycloidal magnet, the incom- 
ing circular polarization finds its handedness naturally 
coupled to the sense of rotation of the magnet moments, 
allowing the population of domains to be determined by 
observing differences in the intensity of magnetic diffrac- 
tion peaks under illumination with LCP or RCP X-rays. 
An Au (222) crystal analyzer system mounted on the 
detector arm of a six-circle diffractometer was used to 
characterise the polarization of the scattered photons. 
The polarization is represented by the Stokes vector 
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FIG. 2: (a) Schematic of experimental scattering geome- 
try used to determine the polarization dependence of A-type 
magnetic satellite reflections (4, ±r, ±1) in TbMnOa. (b) In- 
cident polarization dependence of these reflections in the FE 
cycloidal phase at 15 K, obtained after annealing the sample 
with E < 0. The polarization of the scattered photons was 
analyzed in the tt' channel (r; — 90°). 



P={P^,P2,P3) [2l|,l2g], whose linear components Pi, P2 
were determined by fitting the dependence of the scat- 
tered X-ray intensity 1(77) — Io{l + Pi cos 2rj + P2 sin 2?^) 
on the angle rj of the polarization analyzer about the 
wavevector k'. The Stokes parameters of the incident 
beam were carefully checked before and after every set of 
scans by using the same polarization analyzer, to control 
the incident circular photon polarization, which is very 
sensitive to the beam position stability. 

The initial objective of our study was to establish 
whether or not our experimental setup provided suffi- 
cient sensitivity to any imbalance in the population of 
magnetic domains produced by the applied E field. Fig- 
ure [21(b) shows the polarization dependence of the scat- 
tering from the four A-type [2^ (4, ±t, ±1) magnetic 
satellite refiections — each satellite forming one arm of 
the so-called star of wavevectors around (4,0,0), see 
Fig. [IKa). The measurements were performed at 15 K 
in the FE phase, after cooling the sample with a pol- 
ing voltage of —700 V. Data were collected with linear 
(tt, black diamonds), RCP (e+, open green circles), and 
LCP (e~, red filled circles) polarized X-rays, while the 
scattered beam was analyzed in the tt' channel (77 = 90°). 
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FIG. 3: Variation with analyzer rotation angle i] of the X- 
ray magnetic scattering in TbMnOa at 15 K, for field cooling 
with either E > (upper panel) or E < (lower panel). 
The dashed lines represent the model described in Ref. [3|, 
whereas the continuous lines are our model calculations, as 
described in the text. 



Inspection of Fig. [Hb) reveals that the intensities of the 
magnetic satellites are very similar when incident linear 
polarization n is selected. For incident circular polariza- 
tions this is not the case: instead the intensities display 
complementary behaviour, depending on the sign of r. 
In the case of equi-populated cycloidal domains the in- 
tensities associated with e+ and should be similar, 
whereas for a single cycloidal domain a large difference 
is expected, as is indeed observed to be the case. 

The standard approach to obtain a microscopic un- 
derstanding of the field-induced magnetic domain state 
would be to collect intensity data for a set of satellites 
(h,k ± T,l) with (hkl) taking on as wide a range of val- 
ues as possible. The disadvantages of this method are 
that large corrections must be applied for effects such 
as absorption, etc., and that our specific scattering ge- 
ometry restricts the accessible (hkl) values. Instead we 
undertook a detailed analysis of the polarization of the 
scattered beam by carefully measuring 1(77) for each of 
the satellites shown in Fig.[2fa), extracting the Stokes pa- 
rameters. Figure [3] summarises our data for field-coolings 
performed with either E > or E < 0. The data were 
obtained by rocking the crystal analyzer at different 77 
angles, and are normalized by monitor and corrected by 
subtraction of the background, measured off peak. 

Turning first to consider the case of incident linear {tt) 
polarization shown in Fig. [31 it is clear that for a partic- 



ular choice of the direction of E, only small differences 
in 7(77) are observed for the various satellites, and that 
reversing the direction of E does not appear to have any 
effect. Clear differences in /(ry) of the various satellites 
are observed for incident circular polarization. In this 
case reversing the direction of E has a profound effect 
on the observed intensities. For a given field direction, 
changing the sign of r at fixed I leads to a switching 
of the handedness of incident X-rays that produces the 
maximum intensity, whereas the handedness associated 
with the maximum satellite intensity is invariant with 
respect to changes in the sign of L Of considerable sig- 
nificance is the fact that reversing the sign of E leads to 
a switching of the dominant handedness. Also apparent 
are clear mirror symmetries between LCP and RCP X- 
rays for the same reflection, and between ±1 satellites for 
the same incident polarization state. Thus at a qualita- 
tive level the data displayed in Fig. [3] reveal an imbalance 
in the population of the two possible magnetic domains 
created by the applied electric field, an imbalance that 
can be reversed by switching the direction of the field. 

To explore the potential richness of information en- 
coded in the data shown in Fig. [3] we performed exten- 
sive modelling of the magnetic scattering. For the initial 
model we used the structure due to Blasco et al. , as 
shown in Fig. [IJa), and the magnetic structure proposed 
by Kenzelmann et al. 3], as shown in Fig. [TTb) . In Fig. [3] 
the results of our calculations of NRXMS [2l| based on 
this model are represented by dashed lines. For the cal- 
culations the exact scattering geometry of the six-circle 
setup was taken into account. The agreement between 
the model and the data is clearly unsatisfactory. In par- 
ticular, the data show a non zero value of P2 (a linear 
polarization oblique with respect to directions fi', tt'), 
that is not predicted by the model. 

A vastly improved description of the data was achieved 
following the realisation that, in contrast to the earlier 
neutron diffraction study, our experiment yielded sensi- 
tivity to the b component of the ordered moment on the 
Tb sites, while being largely insensitive to the a com- 
ponent. From the fits we are able to fix the phase angle 
between the magnetic modulations of the Tb and the Mn 
atoms. Consider one Mn atom, and the subsequent Tb 
atom moving along ±c (sign determined accordingly to 
the direction of the electric polarization) , then the phase 
angle between their magnetic modulations, evaluated at 
the same coordinate in space, is found to be 1.0 ± 0.1 
TT (with reference to the b components). The phase shift 
between the two orbits of Tb atoms (at z — 1/A and 
z = 3/4 of the unit cell) is determined to be 1.0 ± 0.2 tt. 

The model results in the following NRXMS amplitude 
for an A- type peak of the form (4 ±r ±1): 
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where the value of a = ±1 selects the sign of r, f3 the 
sign of / = ±1 and 7 = ±1 identifies the two cycloidal do- 
mains 1 and 2, respectively (Fig. [T{c)). The vectors A" 
and B contain the dependence on the polarization of the 
incident and diffracted X-rays [2l[. The spin S'i(K) and 
the orbital Li(K) magnetic components contain the form 
factors for Mn^+ and Th^+ at a given scattering vector 
K. We suppose also that the orbital contribution of Mn 
is quenched, and that of Tb equal to its spin moment 
(L'^VS'^'' = 1) as follows from Hund's rules for a ^Fg 
electronic configuration. Finally A^^ describes the frac- 
tional atomic coordinate along the a axis of the Wyckoff 
position (4c), for the Tb atoms [l8|. 

A simplified expression for the NRXMS amplitude in 
the a' and vr' channels, in which we set the scattering 
angle 29 = 90°, A^^ = and place the b axis perpen- 
dicular to the scattering plane, illustrates the sensitivity 
in our experiment to the different components of the Tb 
and Mn magnetic moments, and the dependence of the 
scattered intensity on the handedness of the X-rays, etc: 



NRXMS experiments when circularly polarized X-rays 
are combined with full polarization analysis of the scat- 
tered beam. Applied to cycloidal systems this technique 
is capable of providing a quantitative, microscopic de- 
scription of the cycloidal domain state. For the specific 
case of multiferroic TbMnOa, with two different types 
of magnetic ion, making it by any standards a challeng- 
ing test case, we have shown surprisingly that this ap- 
proach allows us to make important refinements to the 
magnetic structure obtained from neutron diffraction. 
This technique could readily be applied to other mul- 
tiferroics, particularly those such as BiFeOa 24, 2^ for 



which open questions remain concerning the magnetic 
structure of bulk samples and its modification in thin 

films HSIil,!!!. 
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Here e selects the handedness of the incident X-rays, and 
etc., are the projected magnetization densities 
on the incident (scattered) wavevector. The dominant 
Mn terms are real in fs-' and imaginary in fjt', giving a 
scattered beam that is mainly circularly polarized. They 
are large (small) if eaj = +1 (— 1). A non zero value of 
P2 arises when the pair fs-' ± /#' have different moduli. 
This depends on S^^ and Lf''. Their sign in /#/ {fa') is 
to be compared to the sign of the largest term S*^^" in 
fa' (/#'). Then it is seen that adds with opposite 
signs in the real and imaginary parts, so that its effect in 
P2 nearly cancels. Instead Lf' produces a net P2 with 
the sign +€Pj. 

The full expression for /(K) was fitted to the data 
with the outcome represented by the solid lines in Fig. [S] 
It is clear that inclusion of a 6 component of the Tb 
moment leads to an excellent description of the data, with 
the best fit obtained with (L^** + 28^'') equal to 1.0 ± 0.3 
liB- Most importantly, our refined model of the magnetic 
structure allows us to obtain an accurate description of 
the domain state in our sample as a function of applied 
electric field, including the absolute sense of rotation of 
the magnetic moments. We conclude that cooling in a 
positive electric field led to a population of 96(3)% of the 
cycloidal domain in which the transverse spiral of the 
Mn atoms is counterclockwise, when moving along -f b 
and looking from -|-a (domain 1 in Fig.[ljc)), while field 
cooling in a negative electric field produced a population 
of 83(2)% for the domain of opposite sense of rotation 
(domain 2 in Fig. [Ijc)). 

Our results establish the benefits of performing 
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